Dynamics on Spatial Networks and the Effect of Distance Coarse graining 
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Very recently, a kind of spatial network constructed with power-law distance distribution and total 
energy constriction is proposed. Moreover, it has been pointed out that such spatial networks have 
the optimal exponents 5 in the power-law distance distribution for the average shortest path, traffic 
dynamics and navigation. Because the distance is estimated approximately in real world, we present 
an distance coarse graining procedure to generate the binary spatial networks in this paper. We find 
that the distance coarse graining procedure will result in the shifting of the optimal exponents 8. 
Interestingly, when the network is large enough, the effect of distance coarse graining can be ignored 
eventually. Additionally, we also study some main dynamic processes including traffic dynamics, 
navigation, synchronization and percolation on this spatial networks with coarse grained distance. 
The results lead us to the enhancement of spatial networks' specifical functions. 
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PACS numbers: 89. 75. He, 89.75.-k, 89.75.Fb 



I. INTRODUCTION 

The research on complex networks has been one of the 
most active fields not only in physics but also in other 
various disciplines of natural and social sciences [lJ-Q. In 
traditional statistical mechanics, interaction mainly ex- 
ists between neighboring elements. By introducing the 
complex topology of the networks, the whole system can 
emerge some new properties such as small- world [2], scale- 
free degree distribution Q and community structures [J]. 
However, the spatial property is of great significance as 
well, which makes the interaction between nodes go be- 
yond the neighboring effect but under the restraint of 
their underlying geographical site. This property mat- 
ters much in lots of empirical networks including neu- 
ral networkjfj], communication networks [61, the electric- 
power grid[7[, tr ansp ortation svstems|8l4lOl] and even so- 
cial networks flh - KLa ] . Generally, the geography informa- 
tion of the nodes and the distance between nodes in these 
networks would determine the characteristics of the net- 
work and play an important role in the dynamics hap- 
pening in the network. 

About the networks embedded in the geographical 
space, many works has been done pI6l - [23 | . The first cate- 
gory is focusing on the spatial distribution of the nodes 
of these empirical spatial networks flfl fl7| . Specifically, 
networks with strong geographical constraints, such as 
power grids or transport networks, are found with frac- 
tal scaling[16j. Besides, others researchers discussed the 
small-world behavior and the scale-free networks in Eu- 
clidean space [l^, [13. For example, when supplement- 
ing long range links whose lengths are distributed ac- 
cording to q(l) oc l~ a to D-dimensional lattices, Sen, 
Banerjee and Biswas conjectured that the two transition 
points from random networks and Regular networks to 



the networks with small-world effect in any dimension 
are: a = D and a = D + 1 respectively [la]. Also, 
Xulvi-Brunet and Sokolov constructed an growing net- 
work model by f| ex kil~, c 



where Uj is the distance be- 
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tween i and j. Numerical simulations have shown that for 
a < 1 the degree distribution is a power-law distribution 
and for a > 1 it is fitted by a stretched exponential |19| . 
In addition, some researchers also introduced some ways 
to model the empirical geographical networks [la. [2fj| - [22j . 
However, few of these former works in spatial networks 
are related to the total cost restraint. In fact, The total 
cost is very important when designing these real spa- 
tial networks. Because the longer a link is, the more it 
will cost. Very recently, some researches took this aspect 
into account [lj, [23|, [25J. In (23[, based on a regular net- 
work and subject to a limited cost C, long range connec- 
tions are added with power-law distance distribution un- 
der the probability density function (PDF) P(r) — ar~ s . 
Some basic topological properties of the network with 
different 8 are studied. It is found that the network has 
the minimum average shortest path when 8 = 2 in one- 
dimensional spatial networks, In addition, the authors in- 
vestigated a classic traffic model on this model networks. 
It is found that 8 — 1.5 is the optinrization value for the 
traffic process on the spatial networks. In [2q], pairs of 
sites ij in 2-dimensional lattices are randomly chosen to 
receive long-range connections with probability Pr(u,v) 
proportional to r~". With the total energy restriction, 
G. Li et cl. found a — 3 is corresponding to the minimum 
average shortest path. Moreover, they claimed that the 
optimal value for navigation is a = 3 in 2-dimensional 
spatial networks and a = 2 in 1-dimensional ones 24 1. 
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In many empirical researches, especially the study 
on traffic networks, the distances are estimate 
approximately [a-tlCl [29| . That is to say scientists in- 
cline to regard a range of distance as a typical value. 
This procedure named distance coarse graining should 
also be studied in the spatial network models. Actu- 
ally, the coarse graining process has been discussed since 



long times ago [26l - l29j . All the related works focus on 
how to reduce the size of the networks while keep other 
properties unchanged such as deg ree distribution cluster 
coefficient, degree correlation 26] , random walks[27] and 
svnchronizabilitv|28j . On the other hand, the probabil- 
ity distribution of these long range-connections is cho- 



sen as P(r) 



in the spatial network. It means 



most of the connections are short while a few connections 
are relatively long. However, each node has only two 
neighbors. When the total cost constraint is chosen as 
a certain large value in the binary network, the network 
can not provide enough short long-range connections. In 
this paper, we use the distance coarse graining to solve 
this problem. Specially, we will study how the distance 
coarse graining affects the topology and dynamical pro- 
cess in the spatial network model. The result shows that 
the 6 for minimum average shortest p ath, o p timal traffic 
processJlS [H| and navigation [f(| [H 151 137| will shift to 
smaller value. And the more we coarse grain the distance, 
the more significantly the optimal 6 will shift. Interest- 
ingly, when the network is large enough, the effect of 
distance coarse graining can be ignored. In other aspect, 
how the dynamic processes perform in spatial networks is 
an interesting problem but hasn't received enough atten- 
tion. Investigating the dynamics on the spatial networks 
can not only lead to enhancement of the function of the 
spatial network but also provide us with a better under- 
standing of it. Here, we study two more main dynamic 
processes and find synchronizability[38l |39| can also be 
optimized by a typical d while there is no optimal 6 for 
percolation[40| in such spatial network model. 



II. GENERATING BINARY SPATIAL 

NETWORKS WITH COARSE GRAINED 

DISTANCE 

The model network is embedded in a fc-dimensional 
regular network. The long range connections is generated 
from a power-law distance distribution. A total cost C 
is introduced to this network model. Every edge has a 
cost c which is linear proportion to its distance r. For 
simplification, the edge connecting node i and j cause a 
cost represented by its length r^ in the model. 

According to many empirical studies, the distance 
obeys the power-law distribution jllHla I35J . |36I |. Spe- 
cially, for Japanese airline networks, even there is an ex- 
ponential decay in domestic flights, the distance distri- 
bution follows power-law when international flights are 
added[35]. For the U.S. intercity passenger air trans- 
portation network, the distribution of the edge distance 
has a p ower-law tail with the exponent 6 = 2.20 ± 
0.19[lf. Additionally, authors in ref.[IMl found 5=1 
for social systems like the mobile phone communication 
networks. Therefore, the probability distribution of these 
long range connections here is chosen as P(r) = ar~ 
(PDF) in the spatial model. In order to coarse grain 
the distance, when adding a long range link to the net- 



work, we divide the distance into many continuous parts 
in which all the distances are considered as one typical 
value. So the spatial network is constructed as following: 

1. N nodes are arranged in a 1-dimensional lattice. 
Every node is connected with its nearest neighbors 
which can keep every node reachable. Additionally, 
between any pair of nodes there is a well defined 
lattice distance. 

2. Set the approximate interval W of coarse graining 
process. For each node, divide the distance between 
it and other nodes into -^ parts (N max is the 
largest distance between any nodes in the initial 
network), the mth part are (m — 1)W + l,(m — 
l)W + 2,...,mW. 

3. A node i is chosen randomly, and a certain dis- 
tance r(2 < r < N ni ax) is generated with probabil- 
ity P(r) = ar~ , where a is determined from the 
normalization condition ^^^ P( r ) = 1- 

4. Find the part that the distance r belongs to, say 
mth part here. Then, one of the nodes in this part 
is picked randomly, for example node j. An edge 
between nodes i and j is created if there exists no 
edge between them yet. 

5. After step 4, a certain cost r,j is generated. Repeat 
step 3 and 4 until the total cost reaches C. 

Obviously, there are two significant features of the spa- 
tial network: the power-law distribution of the long range 
connections in the network and the restriction on total 
energy. Firstly, we should determine how to choose an 
appropriate W under different total cost C . In the spatial 
network, the probability distribution of these long range- 
connections is chosen as P(r) = ar~ s . It means most of 
the connections are short while a few connections are rel- 
atively long. However, each node has only two neighbors. 
When the total cost reaches a certain value in the binary 
network, the short range part of the network would be- 
come nearly full connected and can not provide further 
short long-range connections. Consequently, binary spa- 
tial network can not be generated directly without the 
distance coarse grained. Clearly, if W is too small, the 
spatial network model still can not provide enough short 
long-range connections. On the contrary, because all the 
nodes in one distance coarse graining part are regarded 
as the same, if W is too large, too many nodes are con- 
sidered in one distance coarse graining part so that the 
power-law distance distribution will be destroyed signifi- 
cantly. 

In order to deduce the formula of appropriate W , we 
consider an extreme condition with 6 equaling to a very 
large positive value. Under this circumstance, all the 
long-range connections are short. So all of them will 
locate in the first distance coarse graining part of each 
node. Consequently, the total energy is C = (2 + 3 + 4 + 
...+W)N. ForC = en, 2+3+4+. .. +W = c where c is the 
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FIG. 1: (Color online) The total energy of two different binary 
spatial networks. The cycle represents the projected binary 
spatial networks in ref.[ll]. The lozenge stands for the binary 
spatial network with coarse grained distance with W — 7. 
The square is the standard total energy which equals C — nx 
c = 1000x30 = 30000. The results are under 100 independent 
realization. 



average energy on each node. After simplification, W 2 
W-2-2c = Q. So we can get W 
W only has to obey W > \^±^=± 
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2 }fioor where the [•] floor 
represents the operation of rounding downward. In this 
paper, we choose W = 4 when c = 10, W = 7 when 
c = 30 and W = 9 when c = 50. 

As mentioned above, the binary spatial network can 
not be generated directly. In ref 23], to get the binary 
spatial network, authors project the weighted spatial net- 
work into unweighted one by imposing all the weight of 
the existing links to 1 , this will lead to losing total energy. 
On the contrary, when the distance is coarse grained, the 
network can provide enough short-term links. Under this 
circumstance, the binary spatial network can be gener- 
ated independent from the weighted one. In Fig.l, we 
compare the total energy of these two different binary 
spatial networks. Clearly, only the network with coarse 
grained distance can satisfy the total energy limit as the 
standard value. 

For the power-law distance distribution, just like we 
discussed above, there is no doubt that it will be de- 
stroyed when the distance is coarse grained. The results 
for distance distribution is reported in Fig. 2. We can see 
that the distance distribution in projected binary spatial 
networks are almost the same as the standard distance 
power-law distribution. On the contrary, in networks 
with coarse grained distance, the distribution is differ- 
ent from the standard one. This is reasonable, because 
the node in each coarse grained part are regarded ho- 
mogeneous. Consequently, the distribution plot is also 
divided into some continuous parts in which the distance 
distribution is normal. However, when we estimate the 
distance with higher scale, which means marking all the 
distance from (m — 1)W + 1 to mW as a specific value 
(to — \)W + 1, the distance distribution become exactly 
the same with the standard power-law, see the green line 



FIG. 2: (Color online) Distance distribution of two different 
spatial networks. The lozenge represents result of the pro- 
jected binary spatial networks in ref. [11]. The square stands 
for the distance distribution under original scale of the binary 
spatial network with coarse grained distance W = 7. The tri- 
angle is the distance distribution under higher scale (See the 
text) of the coarse grained network. The cycle is the standard 
distance power-law distribution with given n = 10000, c = 30 
and 5 = 3. The results are under 10 independent realization. 
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FIG. 3: (Color online) The result for how the optimal 5 
changes as the distance coarse graining procedure, (a) is the 
average shortest path length for different 8 under different 
coarse graining scale, here N = 1000. (b) is how the S for 
minimum average shortest path changes when the network 
getting larger. The results are under 100 independent real- 
ization. 



in Fig.2. 

One of the most important results in the former works 
about spatial networks is that the minimum shortest path 
happens at S = 2 regardless of the total energy |23l l25j . 
Here, we investigate how the optimal 5 changes as the dis- 
tance coarse graining procedure. Fig. 3 shows the result 
for three different level of distance graining. Obviously, 
the larger W we choose to coarse grain the distance, the 
more severely the optimal point shifts to a lower value. 
However, as the size of the network is getting bigger, the 
shifting effect of the distance coarse graining becomes less 
significant. That is to say, notwithstanding the coarse 
graining in distance, the 5 for the minimum shortest path 
still equals to 2 in large networks. 



III. THE DYNAMICS ON SPATIAL 

NETWORKS AND THE EFFECT OF DISTANCE 

COARSE GRAINING 

Dynamics on spatial networks is an interesting topic, 
for example, it can help us obtain the principle to design 
the optimal transportation networks [2 5]. In this section, 
some main dynamics process will be studied on the spa- 
tial networks. Furthermore, most of the time, people 
tend to coarse grain the distance when designing the real 
networks, especially these transport networks. So un- 
derstanding how the distance coarse graining affects the 
function of the network can be not only of great interest 
but also useful. So we will also discuss the effect of dis- 
tance coarse graining on the function of spatial networks. 
Firstly, based on the results in the former works |23l[25|. 
we know that the spatial properties of the network will 
result in optimal 6 for navigation and traffic process re- 
spectively. So how the optimal <5 shifts with the distance 
coarse graining procedure will be investigated. Secondly, 
we will study the synchronizability in this so-called bi- 
nary spatial network with distance coarse grained. Fi- 
nally, the percolation performance in this binary spatial 
network will be studied as well. 



A. The Effect on Traffic Process 

In traffic process [30J, All the nodes embedded on the 
spatial network are treated as both hosts and routers. 
Every node can deliver at most D packets one step to- 
ward their destinations. At each time step, there are 
R packets generated homogeneously on the nodes in the 
system. The packets are delivered from their own origin 
nodes to destination nodes by special routing strategy. 
There is a critical value R c which can best reflect the 
maximum capability of a system handling its traffic. In 
particular, for R < R c , the numbers of created and deliv- 
ered packets are balanced, leading to a steady free traffic 
flow. For R > R c , traffic congestion occurs as the num- 
ber of accumulated packets increases with time, simply 
because the capacities of the nodes for delivering packets 
are limited. 

In fact, the whole traffic dynamics can be represented 
by analyzing the largest betweenness of the network |31j. 
The betweenness of a node is the number of shortest 
path passing through this node. Note that with the in- 
creasing of parameter R (number of packets generated 
in every step), the system undergoes a continuous phase 
transition to a congested phase. Below the critical value 
R c , there is no accumulation at any node in the net- 
work and the number of packets that arrive at node i is 
Rgi/N(N — 1) on average. Therefore, a particular node 
will collapse when Rgi/N(N — 1) > Di, where gi is the 
betweenness coefficient and Di is the transferring capac- 
ity of node i. Therefore, congestion occurs at the node 
with the largest betweenness. Thus R c can be estimated 
as R c = DiN(N - l)/g ma x, where g max is the largest 
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FIG. 4: (Color online) The effect of distance coarse graining 
on the function of the binary spatial networks, (a) is the crit- 
ical value Rc in traffic process for different 8 under different 
coarse graining scale, where N = 1000. (b) shows how the 
optimal 8 for traffic process changes when the network gets 
larger, (c) is the navigation steps for different 8 under differ- 
ent coarse graining scale, where N — 1000. (d) shows how 
the optimal 8 for navigation changes when the network gets 
larger. The results are under 50 independent realization. 



betweenness coefficient of the network. 

Here, we study the traffic dynamics in the spatial net- 
work with coarse grained distance. The results are given 
in Fig. 4(a) and (b). In Fig. 4(a), it is quite obvious that 
there exists an optimal 5 for Rc, which means the trans- 
port capacity reaches its maximum in this spatial net- 
works. From Fig.4.(b), we can clearly see that this op- 
timal S gets closer to 1.5 gradually as the size of the 
network becomes larger. Actually, even losing the total 
energy, the projected spatial network also hasjdie opti- 
mal 6 for traffic process equalling to 1.5 in ref.[23]. This 
implies that the spatial property play an dominative part 
in the traffic dynamics. 



B. The Effect on Navigation 

Another aspect we are going to investigate here is nav- 
igation. In fact, the navigation process in networks is 
based on local information, which is different from the 
shortest path with global information. Hence, the navi- 
gation reflects another ability of the networks. Accord- 
ing to ref. 25j, we choose the navigation strategy as the 
greedy algorithm[37] in this paper. In the former works, 
Kleinberg found that a — 2 for Pr(u, v) oc r~" is the op- 



timal value in the navigation with the greedy algorithm 
in 2-dimensional spatial networks without total energy 
limit [371] . When adding the energy restriction, G. Li et 
cl. found that the optimal value is a = 3 in 2-dimensional 
spatial networks and a = 2 in 1-dimensional ones[25(. 
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FIG. 5: (Color online) The effect of distance coarse graining 
on the synchronizability of the binary spatial networks, (a) is 
the synchronizability R for different 8 under different coarse 
graining scale, where N = 1000. (b) shows how the optimal 
5 for synchronizability changes when the network gets larger. 
The results are under 50 independent realization. 



Here, though we choose the PDF of the distance distri- 
bution as P(r) = ar , S equals to a in one-dimensional 
space J2J]. What interest us most is that how this opti- 
mal value performs when the distance is coarse grained. 
In Fig. 4 (c) and (d), the results are given. The optimal 
5 shifts to a smaller value as the coarse grained interval 
W gets larger. However, when the size of the networks 
is large enough, the optimal 6 will come back to 2 even 
if the distance gets coarse grained, as show in Fig. 4(d). 



C. The Effect on Synchronizability 

Furthermore, we will study the synchronizability in 
this so-called binary spatial network with distance coarse 
grained. The synchronization is a universal phenomenon 
emerged by a population of dynamically interacting 
units. It plays an important role from physics to biol- 
ogy and has attracted much attention for hundreds of 
years. 

In the former works, the analysis of Master Stability 
Function (MSF) allows us to use the eigenratio R = 
Ajv/A2 of the Laplacian matrix to represent the syn- 
chronizability of a network [39]. Hence, we can calculate 
the synchronizability R under different S. In Fig. 5, syn- 
chronizability of the spatial network is enhanced in some 
spccifical exponent 5. Likewise, there is also an optimal 
6 for synchronizability which is 1.5 approximately when 
the network is large enough. 



D. The Effect on Percolation 

We have known that there are two different ways to ob- 
tain the binary spatial networks. The first is to project 
the weighted spatial networks to an binary one, the sec- 
ond is to generate the binary spatial networks directly 
by distance coarse graining procedure. To begin with, 
we compared the percolation performance of the these 
two kinds of networks. In percolation, we consider what 
happens with a network if a random fraction 1 — p of its 
edges is removed. In this bond percolation problem, the 





FIG. 6: (Color online) The results for the percolation perfor- 
mance of (a)projected binary spatial networks and (b) coarse 
grained binary spatial networks. Here, the network size is 
N = 500. The results are under 50 independent realization. 



giant connected component plays the role of the percola- 
tion cluster, which may be destroyed by decreasing p[40(. 
Obviously, different 5 in spatial networks will result in 
different critical parameter p c . The smaller the critical 
parameter p c is, the better the networks perform in per- 
colation. The results for the percolation performance of 
the these two kinds of binary spatial networks are shown 
in Fig. 6. 

From Fig. 6, we can see that the projected binary spa- 
tial networks have an optimal S for percolation while 
coarse grained binary spatial networks do not. Actually, 
when projecting the weighted spatial networks, lots of 
energy gets lost. So in the projected binary spatial, the 
number of total links become smaller as S gets larger 23] . 
This is why the projected binary spatial networks have 
an optimal S for percolation. Actually, when designing 
the binary spatial network, it is supposed to constrained 
by the total cost. In the percolation, we show that adding 
the total energy constraint by the distance coarse grain- 
ing, the spatial network may perform entirely different in 
some functions. It indicates the distance coarse graining 
is necessary when analyzing some dynamics on spatial 
networks. 



IV. CONCLUSION 

The complex networks has been a hot topic in science 
for more than ten years. Works in this field are based 
on the topology of the networks. So far, many empirical 
works claim that the distance distribution of real net- 
works obeys power-law distribution. In theoretical mod- 
eling aspect, researches begin to pay attention to these 
spatial networks. Very recently, these spatial networks 
are constructed with power-law distance distribution and 
under total energy restriction. This kind of spatial net- 
works can reflect the trade off property of real networks 
between the efficiency (distance power-law distribution) 
and the total cost. 

Studying the dynamic on spatial networks is of great 
significance, which can lead to the enhancement of net- 
works' specifical function by choosing the proper expo- 
nent S. In previous works, authors found that there exist 
stable optimal power-law indexes S for minimum shortest 



path, traffic process and navigation. With these under- 
standings, people may obtain the principle to design the 
effective transport networks. In this paper, we study the 
percolation and synchronizability in such binary spatial 
networks with coarse grained distance, we find that syn- 
chronizability can also be optimized by a typical 5 while 
no optimal S exists for percolation in such spatial network 
model. 

In most case of real lives and empirical researches, the 
distances are estimate approximately. In other words, 
people incline to regard a range of distance as a typical 
distance. How this coarse graining procedure affects the 
optimal 5 is also studied in this paper. Our results show 
that the distance coarse graining procedure will make 
the optimal exponent S in power-law distance distribu- 
tion shift to smaller values for all of average shortest 
path, traffic process and navigation. Interestingly, when 
the network is large enough, the effect of distance coarse 
graining can be ignored. As the real networks, say the 



transport networks, is usual of relatively large size, the 
result indicates that the optimal index 5 still works in 
designing principles for the optimal real transport net- 
works. These results above indicate that the distance 
coarse graining can be used as a universal way to generate 
the binary spatial model with its total cost constraint sat- 
isfied and its power-law distance distribution preserved 
effectively. Moreover, investigation of functions of spa- 
tial related networks with coarse grained distance can be 
an interesting extension. 
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